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TeV blazars offer an exciting prospect for discovering cosmological magnetic fields and for probing 
high energy processes, including CP violation, in the early universe. We propose a method for 
reconstructing both the non-helical and the helical magnetic field correlators using observations of 
cascade photons from TeV blazars. 



^ ! I. INTRODUCTION 

o : 

OA ■ Large-scale magnetic fields with micro Gauss strength have been observed in galaxies [1, 2] and clusters of galaxies [3, 
4]. Although it is assumed that such magnetic fields are the result of amplification from weak seed magnetic fields, the 
origin of the seed magnetic field has not yet been understood. There are two classes of models for the seed generation 
^! ' mechanism: astrophysical and cosmological models. In astrophysical models, the seeds are associated with nonlinear 
structures, and are generated during structure formation through the Biermann battery effect [5], or the Weibel 
instability [6]. On the other hand, seeds in the cosmological models are produced in the early universe and exist as 
cxtragalactic magnetic fields. This class of models includes the generation mechanisms in an inflationary epoch, during 
cosmological phase transitions, and at the epoch of recombination (for recent reviews, see [7-9]). 

Recent gamma ray observations suggest the existence of cosmological magnetic fields stronger than ~ 10~ 16 Gauss 
in the voids [10-12]. Although more detailed work will be required to establish this lower limit [13-15], the suggestion 
of such magnetic fields in large-scale structure voids, ~ 100 Mpc away from non-linear structures, provides a strong 
argument in favor of cosmological models. Further support in favor of an early universe origin can be obtained if the 
q ■ cosmological galactic magnetic fields are "helical" , thus indicating a process of magneto-genesis that fundamentally 
$H , violates invariance under charge conjugation plus parity reflection (CP) [16]. 
ryj ■ Magnetic helicity density is defined as 
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where A is the vector potential of magnetic fields B with B = V x A. Magnetic helicity is odd under CP transforma- 
tions as A and B are odd under C, while A is even but B is odd under P. Non-zero magnetic helicity is predicted in 
scenarios in which cosmic baryogenesis and magneto-genesis occur concurrently during a cosmological phase transi- 
, tion [16-20]. Then the magnetic helicity density is related to the cosmic baryon number density and the CP violation 
responsible for the excess of matter over antimatter also provides helicity to the magnetic field. Other scenarios that 
\ can generate helical magnetic fields that are not tied to baryogenesis have been studied in Refs. [21-23]. Helicity is 
£f~) • also an important factor in the evolution of magnetic fields because it helps to transfer magnetic field energy from 
small to large length scales, a process called an "inverse cascade" . Due to the inverse cascade, helical magnetic fields 
can grow to astrophysically relevant scales at the present epoch, even though the initial scale of the magnetic field 
extends only up to the much smaller cosmological horizon scale at the time of the phase transition. 

The two point correlator of a stochastic, homogeneous, isotropic magnetic field contains two independent functions. 
In physical space, these are the normal and helical correlation functions (Mjv(r) and Mji(r)); in momentum space, 
these are the symmetric and antisymmetric power spectra (S(k) and A(k)). There are several tools to measure the 
normal correlator (or the symmetric power spectrum) of cosmological magnetic fields but measuring helicity directly 
is more challenging. The detection of helicity of cosmological magnetic fields has been a primary motivation for this 
work, as few other direct schemes to detect helicity have been proposed [24]. 

A measurement of the normal correlator itself, together with some theoretical input, provides indirect access to the 
helicity of the magnetic field. This is because the presence of helicity influences the evolution of the normal correlator, 
and the exponents characterizing the normal correlator can give us some 1 information about the magnetic helicity (for 
the evolution in the cosmological context, see Ref.[25-30]). The shape of the normal correlator, in particular the 
existence of a peak in the distribution can also inform us about the epoch of magneto-genesis. If magnetic fields are 
produced during inflation, the symmetric power spectrum can be expected to be scale invariant with significant power 
on the present day horizon scale. 

Direct measures of the helical power spectrum arc more difficult. For astrophysical magnetic helicity, it has been 
suggested that the correlation between Faraday rotation measurement and the polarization degree of radio synchrotron 
emission [31-33] can be used. For cosmological magnetic helicity, Rcf. [24] has shown that correlations in the arrival 
momenta of very high energy cosmic rays are sensitive to the intervening magnetic helicity provided the cosmic ray 
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source locations are known. Cosmic Microwave Background (CMB) anisotropy observations may also permit a measure 
of the hclicity through non- vanishing cross-correlation between the temperature and B-mode polarization anisotropics 
and between the E-mode and B-modc polarization anisotropics [34-36]. 

In the present paper, we propose a scheme to measure both the normal and the helical correlation functions of 
cosmological magnetic fields by using TeV blazar observations. Gamma rays with energy greater than ~ 1 TeV can 
scatter with ambient "extragalactic background light" (EBL) photons to pair produce electrons and positrons. The 
generated electrons and positrons create a secondary cascade of GeV gamma rays through the Inverse Compton (IC) 
scattering of CMB photons. As the electron and positron trajectories are bent due to the Lorentz force by a magnetic 
field, the GeV photon cascade carries information about the structure of the cosmological magnetic field. We show that 
cross-correlations between the arrival directions of the secondary cascade gamma rays at different energies are related 
to the correlation function of the extragalactic magnetic fields. If we imagine the cascade photon arrival direction to 
be a vector in the plane of observation, the inner product of vectors at different energies gives the normal part of the 
magnetic field correlation, and the outer product gives the helical correlator. 

In Sec. II we describe the basic geometry of the process and evaluate the arrival direction of cascade photons as a 
function of the intervening magnetic field. In Sec. Ill we evaluate correlators of the arrival directions of cascade photons 
and relate them to the magnetic field correlation functions. Our results are placed in the context of observations in 
Sec. IV where we discuss estimators for the theoretical correlation functions we have found in Sec. III. Our analysis 
uses many simplifying assumptions that we discuss together with conclusions in Sec. V. Throughout this paper, we 
use natural units: h = 1 = c. 



II. GEOMETRICAL SETUP AND DEFLECTION ANGLE 

TeV gamma rays from distant sources at redshift z s cannot propagate freely over cosmological distances because, 
at such energies, interaction with the EBL can produce electrons and positrons. The mean free path of a gamma ray 
with energy i?TeV is given by 

ZW^Tcv) ~ 80^^ M P c (J^) \ (2) 

where k is a numerical factor which accounts for the model uncertainties of EBL. Here we take k ~ 1 [37]. Electrons 
and positrons generated by the TeV gamma ray lose energy by the production of a secondary gamma ray cascade 
through the IC scattering of CMB photons, and have a mean free path, 

^ e ^30kpc(l + z e )- 4 (^ I ^ v ^ , (3) 

where z e is the typical redshift at which TeV gamma rays create pairs, and E e ~ i?T e v/2 is the electron energy. The 
up-scattered CMB photon has energy, 

JS, =|d + *>-W. (£)'- 88 GOV (J^) 2 , (4) 

where ecmb = 6 x 10~ 4 (1 + z e ) eV is the typical energy of CMB photons. Therefore we observe cascade gamma rays 
with energy ~ E 1 , in addition to TeV gamma rays with energy ~ i/rev- 

It is believed that TeV gamma rays are beamed from distant blazars in a narrow jet of opening angle 8j ~ 5° [r/10]" 1 
where T is the Lorentz factor of the gamma ray emitting plasma. The highest energy photons are observed if we lie 
within the opening angle of the jet. Denoting the angle between the source direction and the jet axis by O , we require 
9 Q < 9j. As emphasized in Ref. [38], the most likely situation is that we are located on the edge of the cone, as depicted 
in Fig. 1. 

Suppose that the blazar is located at x s = D s n s in the observer frame. We will assume for simplicity that the 
redshift of the source is less than one and take 1 + z a ~ 1. If there are no magnetic fields in the IGM, we will observe 
cascade gamma rays due to TeV gamma rays also in the direction n s . However if cosmological magnetic fields are 
present, the cascade gamma rays arrive from different directions. We now evaluate the arrival direction assuming a 
stochastic, homogeneous, and isotropic magnetic field. 

Consider an observed cascade gamma ray that resulted from a TeV gamma ray with energy -&rev that was emitted 
at t = at an angle 6 e from the line of sight as depicted in Fig. 1. The TeV gamma ray produces an electron at 
the position Xi at time t = ij, where tj = I?TeV- The produced electron has momentum Pi whose amplitude Pi 
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FIG. 1. The blazar on the left beams TeV photons within a jet of opening angle 0j. The observer is most likely located at 
the edge of the jet, not on the axis. TeV photons pair produce after propagating a distance DtcV- The pairs are bent by 
ambient magnetic fields and up-scatter CMB photons that propagate a distance D 1 to the observer. The emission angle 8 e , 
the observation direction 0, the distance to the source D s , and the pair creation and IC scattering event positions, (xu, Sxi), 
(xLf,5xf) are also shown. 



corresponds to the energy E e w i?Tcv/2. Since the opening angle of the electron-positron pair is very small, of order 
m e /E e ~ 10 -7 , the direction Pi is the same as the direction of the initial TeV gamma ray. 
The momentum of the electron changes on propagation due to the Lorentz force, 

P(t) =P l + qJ^ dt' v(t') x B(x(t')) (5) 

where q = ±e is the electron/positron charge, x(t) and v(t) are the position and velocity of the electron (or positron), 
namely v(t) = x(t) where the overdot denotes differentiation with respect to time. (For convenience, from now on we 
shall refer to the charged particle as being the electron.) 

We will now decompose all vectors in components parallel and perpendicular to the source direction. For example, 
the momentum and the position of the electron at time t is decomposed as 

P{t)=P L (t) + Sp(t), x(t) = x L (t) + 6x(t), (6) 

where the subscript L means the component parallel to the source direction (line-of-sight for TeV source). Therefore, 
the vector Sp(t) and Sx(t) are the deviations induced by the magnetic field. 
In terms of the decomposed components, Eq. (5) can be written as 

P L (t) + 5p(t) = P Li + S Pi +qj dt' [v L (t') + 5v(t')} x B(x(t')), (7) 

where Pu and Spi are the momentum components at time t — ti. Note that at this stage, instead of replacing x by 
Si in the argument of B, we perform the integration along the actual path x(t). This is important if the magnetic 
fields have significant power on small scales, i.e., a blue spectrum. 

The bending angle of the electron is estimated as 6 = D e /R L ~ 1.2 x 10" 3 [B/10" 16 G][£; Tc v/10 TeV]" 2 where 
Rl = E e /qB is the Larmor radius. Here we have assumed a magnetic field coherence scale larger than D e ; otherwise 
the electron trajectory would be diffusive yielding a smaller estimate for 6. Then the maximum deviation from the 
source direction is Sxi ~ 90 kpc(l — DtcV / D S )[B /10~ 16 G][i?Tcv/10 TcV]~ 3 . Since the bending angle is small, we can 
treat 5p, Sx and Sv as perturbations. To linear order in the magnetic field strength, Eq. (7) becomes 

Sp(t) = S P i +qj dt' v L {t') x B(x(t')). (8) 

The electron energy E e is constant during this process, since a magnetic field does no work. Dividing Eq. (8) by 
E e , we obtain the velocity, 



Sv(t) = 5vi + J- J dt' v L (t') x B(x{t')), 



(9) 
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and another integration gives the trajectory, 

Sx[t) - Saa = Svi[t - U) + 4r dt " dt' v L {t') x B(x(t r j). (10) 

&e J ti Jti 

At any point in the electron's trajectory there is a probability ol an IC scattering event with a CMB photon. 
We simplify the present analysis by assuming that the electron travels a fixed distance D e (E e ), and at final time tf 
up-scatters a CMB photon. Since the electron is ultra-relativistic with Lorentz boost factor E e /m e ~ 10 7 , the cascade 
gamma ray propagates along the momentum of the electron. Therefore the arrival direction of the cascade gamma 
ray is the same as the direction of the position Xf = x(tf). 

We now define the vector in the observation plane by 

„ Sxi — Sx f . . 

where 5xj is Sx(t) at t = tf. The magnitude, |0|, is the observed angle shown in Fig. 1, and the direction in the 
observation plane corresponds to the azimuthal direction of the original TeV gamma ray. 

We will now relate to the magnetic field using Eq. (10). However Svi is still unknown in Eq. (10). To evaluate 
it, we first need to express the emission angle 6 e in terms of the observed angle 0. Applying the trigonometric law of 
sines to the triangle formed by the source, observer and electron position - recall that D e <C D s - we obtain 

a ~ D i e « Ds-d™ ^ (12) 



-DtcV £>TeV 

In terms of 0, the vector Svi can now be written as 

5v t = v e 6 e y = v e — DtcV 0, (13) 

i^TcV 

where v e is the magnitude of the electron velocity, and y is the unit vector perpendicular to the line-of-sight (see 
Fig. 1). 

Now Eqs. (10), (11) and (13) give, 



e(£r) = J'' dt " [ dt ' vlV) x B(x(t% 



(14) 



where we can rewrite D^ c y, D e , and E e as functions of E, 
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£tcv(£tcv) ~ 80 Mpc ( ) , (15) 



-1/2 



°e - 30 kpc ( ) , (16) 



E-, 



1/2 



Ee - 10 TeV ( ) . 



With these relations, the magnitude of 0(£? 7 ) is roughly estimated as 

^ - " — - >< ^ (iWksi) (lOb^ev) (ib^y 



(18) 



III. CORRELATORS AND PREDICTIONS 

The vector &(E y ) describes the position of the observed cascade gamma ray in the observation plane but its 
magnitude depends on geometrical factors such as the distance to the source (Eq. (14)). We define a rescaled vector 
to remove such dependence, 

Q(£ 7 ) = ^ Ds 0(£ 7 ). (19) 
qD Tc yD e 
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Now we are interested in two types of correlators of the Q(E 1 ) vectors 

F(E U E 2 ) = {Q{E 1 )-Q{E 2 )), 



(20) 



G{E U E 2 ) = (Q(E 1 ) x Q{E 2 )-x) 



(21) 



where E\ and E 2 are two energies of the observed cascade gamma rays, and we recall that we have set up our 
coordinate system so that the x— axis is along the line of sight, i.e., x is normal to the observation plane. 

The ensemble value (Qi(Ei)Qi' (E 2 )) , where i and i' denote components in the observation plane, can be written 
by using the magnetic field correlation function, Eq. (14), 



(Q l (E 1 )Q l ,(E 2 )) = eyje^i^iK (^) 



dt[ f*' 1 dt x 



u D e i Jtn D e i 



dt'o 



dU 



-L / (B\x(t 1 ,E 1 ))B l (x(t 2 ,E 2 ))), (22) 

U e.1 Jti U e2 



where D e \ and D e2 are D e for the cascade gamma ray with energy Ei and E 2 respectively. Note that the velocity 
vl = (1, 0, 0) and tf — t\ = D e . 

The correlation function of a stochastic, homogeneous, and isotropic magnetic field is given by [39] 



(B i (x + r)B j (x))=M N {r) 



+ M L (r)^ + M H (r)e ljl r l , 



(23) 



where M^v(r), M[ J (r) 1 and Mh{t) are the correlation functions for the normal, longitudinal, and helical parts of the 
magnetic fields. Due to the homogeneity and isotropy of the magnetic fields, these correlations depend only on the 
separation distance r = \r\. The divergence- less condition gives 



M N {r) 



1 d 

2r dr 



(r 2 M L (r)). 



(24) 



The stochastic fields are often described by power spectra in Fourier space. The magnetic field correlation function 
in Fourier space is 



(B:(k)B J (k / ))^(2nfS 3 (k-k / ) 



(25) 



where S(k) and A(k) are the symmetric and antisymmetric (helical) parts of the magnetic field power spectrum. The 
functions S(k) and A(k) are related to the correlation functions Mjv(r), M^(r), and Mn{r) as in [39]. 

Therefore the ensemble average, Eq. (22), can also be decomposed into three parts: the normal, longitudinal, and 
helical parts, 



(Qi(Ei)Qi'(E 2 )) — Cnw + Cm' + Chw 



(26) 



with each of the terms given by four integrations as in Eq. (22). The separation scale, r, appearing in the correlators 
AlN(r), Ml(t) and Mh{t) is the magnitude of the separation vector 



r(t 1 ,t 2 ,E u E 2 ) = x(h,Ei) - x{t 2 ,E 2 ). 
Since cosmological magnetic fields are weak, this can be approximated as 

r(t l7 t 2 ,Ei,E 2 ) ps XL(ti,Ei) - x L (t 2 ,E 2 ). 

Therefore, 

r(t 1 ,t 2 ,E 1 ,E 2 ) = (Dtcv(Ei) - D Tc v(E 2 ) + (h-t u ) - (t 2 -t 2i ))x. 
For ti = tu and t 2 = t 2 {, the separation scale becomes 

r{t lt ,t 2i ,E 1 ,E 2 ) = D TcV (£ x ) - £tcV [E 2 ) + D e {E x ) - D e [E 2 ) . 



(27) 
(28) 
(29) 
(30) 



In the case with i = i' ^ x, Clu and Chu vanish, because only the x components of and r are non- vanishing. 
However the normal correlator does not vanish and Cm a is given by 



Cnu(Ei,E 2 ) 



d^f'dti f dt' 2 f 2 dt 2 Mjfirih^^Ea)). 

tli Jt u Jt 2i Jti 



(31) 
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FIG. 2. Events at two different energies sample the magnetic field in regions of size D e ~ 30 kpc (solid lines at the vertices of 
the triangles). The regions themselves are separated by distance r which can be ~ 100 Mpc depending on the energy difference 
of the two events. Energy resolution of the detector translates into a lower limit on the separation at which the magnetic field 
correlations can be probed. 



In the case with i ^ i' , i ^ x and i' ^ x, while Cnw and Clw vanish, Chw has non-zero value, 

CW = ew» [ 11 dt[ f 1 dh f * dt' 2 f 2 dt 2 M H (r(t 1 ,t 2 ,E 1 ,E 2 ))(r(t 1 ,t 2 ,E 1 ,E 2 ) ■ &), 

Jtii Jtu Jt 2 i Jt 2i 

Therefore F and G can be written directly in terms of the magnetic field correlation functions, 



* lf dt[ r*' 1 dh f t2t dt'i 



dh 



(32) 



(33) 



G(E 1 ,E 2 )=2 



dt[ 



D, 



dh 



el Jt u ^el Jt 2 



—^M H (r{t 1: t 2 ,E 1 ,E 2 ))(r(t 1 ,t 2 ,E 1 ,E 2 )-x). (34) 



It is worth pointing out that F(E y ,E 1 ) ^ but G(E 7 ,E 7 ) — 0. 

We have already seen that -DtcV ^ D e independently of £L, and so the separation scale in the magnetic field 
correlation function can be approximated as, 



r(ti, t 2 , E u E 2 ) w D TeY (Ex) - D TcV (E 2 ). 
Then recalling tf — U = D e , we can do the integrations in Eqs. (33) and (34) to get 

F{E X ,E 2 ) a ht N (\r 12 \), 



(35) 



(36) 



G(E 1 ,E 2 ) a -M H (\r 12 \)r 12 , 



(37) 



where r\ 2 is r\ 2 — Dtcv(Ei) — DTev(E 2 ). This shows that we can obtain both the non-helical and the helical magnetic 
field correlation functions through the correlators F(E 1 ,E 2 ) and G(E 1 ,E 2 ) constructed from the arrival information 
of cascade photons. This is the main result of this paper. 

Note, as shown in Fig. 2, that the obtained correlator is on the scale r X2 which can be much larger than the scale 
D e . This is because the magnetic field gets correlated at the spatial points where pair production occurs and these 
can be separated by hundreds of Mpc depending on the values chosen for E\ and E 2 . It remains an interesting open 
question if further refinements of the above method can lead to correlators on scales smaller than D e ~ 30 kpc. 



IV. OBSERVATIONS AND ESTIMATORS 



In the previous section, we have considered correlators in the arrival directions of cascade photons and expressed 
them in terms of correlators of the intervening magnetic field. In this section we consider the problem of determining 



7 



magnetic field correlators from the observers point of view: given cascade arrival direction data, what quantity should 
be calculated that corresponds to the magnetic field correlator? The calculated quantity will at best be an estimator 
for the magnetic field correlator because observations are limited in number, while the magnetic field correlator in 
Eq. (23) is over an infinite ensemble of realizations. 

In a more realistic setting - see simulations in Ref. [40] - observed gamma rays with some energy E are expected 
to be scattered around a typical observation direction on the observational plane. The scattering is due to differences 
in the magnetic field along the trajectories of different electrons and also the stochasticity of the EBL and CMB 
photons. We have not considered the latter fluctuations in this paper but the variation in the magnetic field can 
be smoothed out by taking the average (&(E)) over all observed gamma rays with the same energy E. Even if the 
magnetic fields have small-scale structure, i.e., a blue spectrum, the averaging procedure should yield the correct 
magnetic field correlator on scales larger than the smoothing scale, \Sx\ ~ 90 kpc. 

There is a second way to smooth out fluctuations in the magnetic field. This is by realizing that the correlation 
functions F{Ex,E 2 ) and G(Ei 1 E 2 ) evaluated in the previous section depend on only one function, r(ti, t 2 , E\, E 2 ). 
Hence there are many choices of E\ and E 2 yielding the same r, and one of the energy variables can be integrated 
out. This corresponds to averaging over all pairs of energies such that the distance r remains fixed (see Fig. 2). Using 
Eq. (35) we see that r(ti,t 2 , E 1 ,E 1 + SE) is independent of E 1 provided 



5E(E yi r) = 27 GeV 



GcV 



10 Mpc 



E- 



GcV 



E- 



(38) 



We can now average Eqs. (36) and (37) over E\ while taking E 2 = E\ + SE : 

/dE 
-^Q{E 1 )-Q{E 1 + 5{E 1 ,r)), 



(39) 



-Af ff (r) 



dE 7 
~AE 



Q(£ 7 ) x Q(E y + S{E^,r)) ■ x, 



(40) 



where AE is the integration range of the observation energy. The vector Q(E) now denotes the average (rescaled) 
direction vector for cascade photons with energy E. 

A third way to perform the ensemble average is to use observations of many TeV blazars because cascade gamma rays 
from different blazars sample magnetic fields along different path. Hence we can obtain the magnetic field correlation 
function using 



M N (r) 



9 N 

-Y 



dE-y 

~Ke 



(41) 



2 N 



dE 1 
~AE 



Q(E 7 ) x Q{E 1 +8{E 1 ,r)) ■ x 



(42) 



where a labels the blazar. 

Finally we would like to remark that the most likely observation of a TeV blazar is when we are positioned at the 
edge of the jet [38], but we are even more likely to be located outside the jet opening angle. Then we will not observe 
the TeV source but will still receive some cascade GeV photons. In this case, if there is reason to suppose that the 
observed photons are indeed from a cascade, we can extend our correlator by replacing Q(E) in Eqs. (41), (42) by 
Q(E) — Q(E*) where E* is the highest energy (and least deviated) cascade photon that is observed. We can still draw 
some conclusions about the magnetic field correlators though the analysis is more involved. 



V. DISCUSSION AND CONCLUSIONS 



Our main result is the connection between cosmological magnetic field correlators and the correlators of cascade 
photons as given in Eqs. (36) and (37) together with the definitions in Eqs. (20) and (21). In the observational context, 
this connection can be written as in Eqs. (41) and (42). These relations suggest that observations of cascade photons 
may be used to directly study the non-helical and helical spectra of cosmological magnetic fields. There are few other 
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direct ways to probe the helicity of cosmological magnetic fields and this is an important feature of the technique we 
have described. 

Our analysis will need further development as it has used many simplifying assumptions. For example, the stochastic 
interactions of the TeV photons with the EBL photons will lead to a probability distribution for -DTeV while we have 
used a fixed value. Similarly there will be probability distributions for the other interactions, and the magnetic held 
will have some structure on small scales. In addition, if the blazar is at high redshift, cosmological expansion will 
become important. These effects can best be studied using Monte Carlo simulations similar to the ones described 
in Refs. [40, 41]. The simulations should be able to confirm if the magnetic field correlators can be recovered from 
realistic data, which hopefully will become available in the near future. 
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